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DYNAMIC  RESPONSE  OF  A  FLUID-LOADED  PLATE 
CONTAINING  PERIODIC  MASSES 


1.  INTRODUCTION 


Plate  theory  has  been  analyzed  for  many  years.  Thin  plate  theory1  is  a  simplified  model 
that  does  not  incorporate  rotary  inertia,  shear  effects,  or  other  plate  dynamics.  In  contrast,  thick 
plate  theory2  usually  incorporates  all  of  the  dynamics  of  the  plate  and  is  normally  used  when  the 
plate  thickness  is  on  the  order  of  a  wavelength  of  propagating  flexural  waves  within  the  plate. 
More  thorough  investigations  have  analyzed  the  propagation  of  waves  (dispersion  curves)  for  the 
plate  without  fluid  loading,3  and  for  the  plate  in  contact  with  a  continuous  fluid  on  one  or  both 
sides.4  6  Recently,  closed-form  solutions  to  the  equations  of  motion  were  derived  for  an  infinite 
thick  plate  coupled  on  one  or  both  sides  with  fluid  loading  as  it  is  excited  with  a  continuous 
forcing  function.7  These  previous  studies'  7  provided  a  thorough  basis  for  understanding  plate 
theory;  however,  they  did  not  take  into  account  the  addition  of  attached  structural  components. 

The  inclusion  of  periodic  mass  or  stiffness  elements  in  structures  has  been  used  for  years  as 

a  method  to  obtain  the  desired  static  and  dynamic  response  of  a  system.  Numerous  investigations 
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of  this  system  have  been  undertaken  to  determine  the  farfield  acoustic  radiation  pattern.  ' 
Throughout  the  past  studies,  a  model  was  loaded  with  a  line  or  point  source  and  the  resulting 
acoustic  field  was  predicted.  The  problem  that  has  largely  been  ignored  is  the  displacement  field 
of  a  thick  plate  containing  periodic  masses  or  stiffeners  when  it  is  excited  by  an  incoming  acoustic 
plane  wave.  In  a  thin  or  thick  plate  with  embedded  masses,  the  displacement  and  stress  fields  are 
subjected  to  different  forces  than  they  would  otherwise  encounter  in  a  homogeneous  plate. 

This  report  derives  the  equations  of  motion  of  a  fluid-loaded  thin  plate  that  contains 
embedded  masses  when  it  is  forced  by  an  acoustic  plane  wave  at  a  specific  wavenumber.  The 
thin  plate  equations  of  motion  are  coupled  to  the  fluid  and  masses  and  then  manipulated  to 
determine  the  displacement  field.  This  method  is  then  expanded  to  Timoshenko-Mindlin 
equations  of  motion,  as  well  as  to  thick  plate  theory  using  finite-element  solutions.  Finally,  the 
thick  plate  equations  of  motion  containing  bottom  attached  masses  are  formulated  and  the 
resultant  problem  is  discussed. 
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2.  PLATE  MODELS  WITHOUT  MASSES 

To  facilitate  a  basic  understanding  of  the  effects  of  periodic  masses,  it  is  first  desirable  to 
list  the  thin  and  thick  plate  theories  without  the  masses  and  graphically  compare  the  solutions  of 
each  plate  model  in  the  wavenumber-frequency  domain.  The  three  plate  models  listed  in  this 
section  have  all  been  previously  derived;  therefore,  only  the  governing  partial  differential 
equations  and  the  final  solutions  are  shown  here.  They  correspond  to  a  thin  plate,  a  thin  plate 
that  incorporates  shear  and  rotary  inertia  effects  (i.e.,  the  Timoshenko-Mindlin  plate),  and  a  thick 
plate.  A  diagram  of  the  modeled  system  is  shown  in  figure  1 .  Note  that  using  the  coordinate 
system  in  the  orientation  shown  in  figure  1  results  in  a  negative  value  of  a\  the  thickness  h  is  a 
positive  value. 

Fluid 


Acoustic  Load  P, 
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2.1.  THIN  PLATE  WITHOUT  MASSES 

The  first  model  investigated  is  that  of  a  thin  plate  with  a  fluid  load  on  top.  The  equation  of 
motion  of  the  thin  plate  is1 1 


d  u  (x,t)  dzu(x,t) 

D - - +  ph - *_ - =  -p(x,t)  , 


dx 


av 


0) 


where  uz(x,t)  is  the  displacement  of  the  plate  in  the  z-direction,  x  is  the  spatial  position,  t  is  time, 
p  is  the  density  of  the  plate,  p(x,t)  is  the  applied  pressure  in  the  fluid  load,  and  D  is  the  bending 
stress  expressed  as 


D  = 


12(1  -o2)  ’ 


(2) 


where  E  is  Young’s  modulus  of  the  plate  and  v  is  Poisson’s  ratio  of  the  plate.  The  acoustic 
pressure  in  the  fluid  is  governed  by  the  wave  equation  and  is  written  as 


^2/?(x,z,Q  <?2/7(x,z,/)  1  g?2p(x,z,/) 

^z2  ^x2  c)  at1 


(3) 


where  p{x,z,t )  is  the  pressure  in  the  fluid  and  cy  is  the  compressional  wavespeed  in  the  fluid.  The 
interface  at  the  fluid  and  solid  surface  of  the  plate  at  z  =  b  satisfies  the  linear  momentum  equation 
and  is  written  as 


Pf 


2 

a  Uz(x,b,t) 


apjxjyp) 

az 


(4) 


where  pf  is  the  density  of  the  fluid.  Assuming  harmonic  response  and  loading  in  space  and 
time,  equations  (1),  (3),  and  (4)  can  be  combined.  The  resulting  expression  is  the  displacement 
in  the  z-direction  divided  by  the  incident  pressure  Pj  and  is  written  as 


U7(kr,co) 


-2 


Dk*  + 


P_L 

iy 


-  ph 


(O 


(5) 


4 


where  kx  is  the  wavenumber  with  respect  to  the  x-axis,  co  is  frequency,  and  /  is  the  acoustic 
wavenumber  written  as 


r= 


(6) 


2.2.  TIMOSHENKO-MINDLIN  PLATE  WITHOUT  MASSES 

The  second  model  investigated  is  that  of  a  thin  plate  that  includes  the  effects  of  rotary 
inertia  and  transverse  shear,  generally  called  a  Timoshenko-Mindlin  plate.  This  plate  is  also  in 
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contact  with  a  fluid  load  on  top.  The  equation  of  motion  for  this  plate  is 


D 


3  x4 


( 


Dp  i  ph 
kG 


3  A  j)4 


12 


3  uz(x,t) 


„2.3\  ^4 


3x23t 2 


pGr 
12  kG 


3  u7(x,t )  3  w,(x,0 


3t 


+  ph- 


3P 


=  -  p(x,t)  + 


f  D  ] 

32  p(x,t) 

(ph2) 

32  p(x,t) 

\xGh; 

3x2 

v  1 2  kG  ) 

3t 2  ’ 

(7) 


where  G  is  the  shear  modulus  of  the  plate  and  /cis  the  shear  coefficient  of  the  plate,  typically 
computed  for  a  rectangular  section  using 


5 

K  = -  . 

6-o 


(8) 


Assuming  harmonic  responses  in  space  and  time,  equations  (3),  (4),  and  (7)  can  be  combined. 
The  resulting  expression  is  the  displacement  in  the  z-direction  divided  by  the  incident  pressure 
and  is  written  as 


U ^(k  ,co) 


-2- 


2D 

K2  + 

(2ph2 ) 

2 

CO 

\  xGh. 

) x 

K 12  kG  j 

pI 

4  (°Pf 
Dk*  + 

Dp 

ph*  'j  2.2 
-y  60  kz.  + 

V*3 

Pfph2 ) 

co4  + 

( Pf  \ 

^  iyKGh 

kG 

12  J 

\2kG 

V 

12iynG 

l  ir  J 

.(9) 


(O 
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2.3.  THICK  PLATE  WITHOUT  MASSES 


The  third  model  investigated  is  that  of  a  thick  plate — this  model  is  sometimes  referred  to  as 
the  exact  or  complete  solution.  The  equation  of  motion  for  a  three-dimensional  medium  is13 

//V2u  +  (T  +  //)VV  •  u  =  p—~-  >  (10) 

dt 2 

where  A  and  / u  are  the  Lame  constants,  •  denotes  a  vector  dot  product;  and  u  is  the  Cartesian 
coordinate  displacement  vector  of  the  plate.  Assuming  harmonic  response  in  space  and  time, 
equations  (3),  (4),  and  (10)  can  be  combined.  The  resulting  expression  is  the  displacement  in  the 
z-direction  divided  by  the  incident  pressure  and  is  written  as7 

£/2(£y,z,ra)  uIs(kx,z,co) 

P,  /jAs(kx,a>)  *  (U) 


where  z  is  the  spatial  location  (in  meters),  and 

Uls  (kx,z,co)  =  /j  {sin(orz)  -  cos (J3h) sin[or(z  +  /?)]} 

+  /2  {sin(/?z)  -  cos(or/?)sin[/?(z  +  /z)]} 

+  /3  sin(«/?)  cos[/?(z  +  h)]  +  /4  sin (fih)  cos[or(z  +  //)],  ( 1 2) 


and 


A5(A:v,ry)  =  f 5  (k  x,  (o)  cos(a  h)  cos(  J3  h)  +  f6(kx,co)cos(ah)sin(/3h) 

+ /7  (kx ,  (o)  sin(«  h )  cos(/?  h)  +  /8  (kx ,  c o )  s'\n(ah)  sin(  J3h)  -f5(kx,o)).  (13) 

The  constants  in  equations  (12)  and  (13)  are 


f\(kx,G>)  =  -8or2/?£2(/?2  -kx)  , 

(14) 

f2(kx,a>)  =  -4cckl(/32-k*)2  , 

(15) 
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f^{kx,co)  =  16 a2pkx  , 

(16) 

f4(kx,o))  =  2a(/32  -k2)3  , 

(17) 

f5(kx,o))  =  -Saj3k2(/32  -k2)2  , 

(18) 

f6(kx,a)  =  ipf(ypY'a(p2  -*2)2(/?2  +i2)2  , 

(19) 

f1{kx,co)  =  *iPf{rpTlci2pk2x{p2  +  k2x)2  , 

(20) 

and 

/8  ( k x ,  ®)  =  (P 2  -  k  2  )4  +  1 6a  2 p  2k  4X  . 

(21) 

In  equations  (12)  through  (21),  a  is  the  modified  wavenumber  associated  with  the  dilatational 
wave  and  is  expressed  as 

a=K-k2x  , 

(22) 

where  kj  is  the  dilatational  wavenumber  and  is  equal  to  cole with  cj  being  the  dilatational 
wavespeed  (m/s);and  /?  is  the  modified  wavenumber  (rad/m)  associated  with  the  shear  wave  and 
is  expressed  as 

(23) 

where  ks  is  the  shear  wavenumber  (rad/m)  and  is  equal  to  co/cs ,  with  cs  being  the  shear 
wavespeed  (m/s).  The  relationships  between  wavespeeds  cj  and  cs  and  the  plate’s  Lame 
constants  X  and  fj  are  determined  by 

1 A  +  2u 

Cdi  P  ■ 

(24) 
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and 


c 


s 


(25) 


Finally,  the  relationships  between  the  Lame  constants  and  Young’s  modulus,  shear  modulus,  and 
Poisson’s  ratio  are  given  by 


A  = 


Eo 

(l  +  u)(l-2u) 


(26) 


and 


/u  =  G  = 


E 

2(1 +  u) 


(27) 


2.4.  COMPARISON  OF  MODELS  WITHOUT  MASSES 

The  models  without  masses  are  now  compared  using  a  numerical  example.  The  properties 
of  the  plate  and  fluid  are  listed  in  table  1 .  The  results  are  displayed  graphically  as  magnitude  of 
normalized  velocity  divided  by  incident  pressure  (a  dimensionless  quantity)  versus  wavenumber. 
This  expression  is  determined  by 


V"(kx,co)  .  Uz(kx,a>) 

- =  '°PfCf - » - ' 


(28) 


Figure  2  is  an  image  of  thin  plate  normalized  velocity  versus  wavenumber  and  frequency  and 
was  generated  using  equations  (5)  and  (28).  Figure  3  is  an  image  of  thin  plate  (with  rotary 
inertia  and  shear  effects)  normalized  velocity  versus  wavenumber  and  frequency  and  was 
generated  using  equations  (9)  and  (28).  Figure  4  is  an  image  of  thick  plate  normalized  velocity 
versus  wavenumber  and  frequency  and  was  generated  using  equations  (1 1)  and  (28).  Figures  5 
through  24  are  constant  frequency  cuts  of  Figures  2,  3,  and  4  overlaid  for  a  frequency  range  of 
500  Hz  to  10  kHz  in  500-Hz  increments.  In  figures  2  through  24,  the  velocity  is  displayed  on  a 
decibel  scale.  Note  that  the  thin  plate  and  thick  plate  solutions  begin  to  diverge  at  about  1 .5  kHz, 
and  the  Timoshenko-Mindlin  plate  and  thick  plate  solutions  begin  to  diverge  at  about  3.5  kHz. 


8 


Wavenumber  (rad/m) 


Table  1.  Parameters  Used  for  Fluid-Loaded  Plate  Without  Mass 


Plate  Properties 

Thickness,  h 

0.04  m 

Young’s  Modulus,  E 

7.0e8  N/m2 

Shear  Modulus,  G 

2.5e8  N/m2 

Poisson’s  Ratio,  v 

0.4  (dimensionless) 

Density,  p 

1 200  kg/m3 

Fluid  Properties 

Compressional  Wavespeed,  Cf 

1500  m/s 

Density,  pp 

1 025  kg/m3 

40 
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Figure  2.  Normalized  Velocity  Versus  Wavenumber  and  Frequency  Without  Masses — 

Modeled  Using  Thin  Plate  Theory 
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Figure  3.  Normalized  Velocity  Versus  Wavenumber  and  Frequency  Without  Masses — 
Modeled  Using  Timoshenko-Mindlin  Plate  Theory 
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Figure  4.  Normalized  Velocity  Versus  Wavenumber  and  Frequency  Without  Masses — 

Modeled  Using  Thick  Plate  Theory 
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ure  5. 


Normalized  Velocity  Versus  Wavenumber  Without  Masses  at  500  Hz 


Figure  6.  Normalized  Velocity  Versus  Wavenumber  Without  Masses  at  1.0  kHz 
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Normalized  Velocity  Versus  Wavenumber  Without  Masses  at  1.5  kHz 


Figure  8.  Normalized  Velocity  Versus  Wavenumber  Without  Masses  at  2.0  kHz 


Magnitude  (dB)  Qg'  Magnitude  (dB) 


ire  9. 


Normalized  Velocity  Versus  Wavenumber  Without  Masses  at  2.5  kHz 


Figure  10.  Normalized  Velocity  Versus  Wavenumber  Without  Masses  at  3.0  kHz 
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ire  11.  Normalized  Velocity  Versus  Wavenumber  Without  Masses  at  3.5  kHz 


Figure  12.  Normalized  Velocity  Versus  Wavenumber  Without  Masses  at  4.0  kHz 
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ire  13.  Normalized  Velocity  Versus  Wavenumber  Without  Masses  at  4.5  kHz 


Figure  14.  Normalized  Velocity  Versus  Wavenumber  Without  Masses  at  5.0  kHz 
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ire  15.  Normalized  Velocity  Versus  Wavenumber  Without  Masses  at  5.5  kHz 
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Figure  16.  Normalized  Velocity  Versus  Wavenumber  Without  Masses  at  6.0  kHz 


Magnitude  (dB)  Magnitude  (dB) 


r re  17.  Normalized  Velocity  Versus  Wavenumber  Without  Masses  at  6.5  kHz 


Figure  18.  Normalized  Velocity  Versus  Wavenumber  Without  Masses  at  7.0  kHz 
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ire  19.  Normalized  Velocity  Versus  Wavenumber  Without  Masses  at  7.5  kHz 


Figure  20.  Normalized  Velocity  Versus  Wavenumber  Without  Masses  at  8.0  kHz 


Figure  21.  Normalized  Velocity  Versus  Wavenumber  Without  Masses  at  8.5  kHz 
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Figure  22.  Normalized  Velocity  Versus  Wavenumber  Without  Masses  at  9.  0  kHz 
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Figure  23.  Normalized  Velocity  Versus  Wavenumber  Without  Masses  at  9.5  kHz 


Figure  24.  Normalized  Velocity  Versus  Wavenumber  Without  Masses  at  10.0  kHz 
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3.  PLATE  MODELS  WITH  PERIODIC  MASSES 


The  three  plate  models  (thin  plate,  Timoshenko-Mindlin  plate,  and  thick  plate)  are  now 
modified  to  account  for  periodic  masses,  which  act  as  applied  forces;  see  the  right-hand  side  of 
equations  (1),  (7),  and  (10).  A  diagram  of  the  modeled  system  is  shown  in  figure  25. 


Fluid 


Acoustic  Load  P, 


Figure  25.  Plate  with  Masses  Subjected  to  a  Fluid  Load 


3.1.  THIN  PLATE  WITH  MASSES 

The  first  model  investigated  is  that  of  a  thin  plate  containing  the  periodic  masses  with  a 
fluid  load  acting  on  the  upper  surface  of  the  plate.  The  equation  of  motion  of  the  thin  plate  is 


d  u  (x,t)  dzu  (x,t) 
D - l-r-  +  ph - - 


dx 


dr 


=  -p(x,t)~  pM(x,t)  , 


(29) 


with 
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(30) 


~  dlu.{x,t) 

pM{x,t)  =  - ~2 - S(x-nL)  , 

dt2 


n=—  oo 


where  M  is  the  mass  per  unit  length  of  the  masses,  L  is  the  spacing  between  adjacent  masses,  and 
8  is  the  Dirac  delta  function.  The  problem  is  transformed  into  the  wavenumber-frequency 
domain  using 


and 


u7(x,t)  =  U  z(kx,(o)exp(ik  Xx)exp(i  cot). 

(31) 

p(x,t)  =  P(kx,co)exp(ikxx)exp(icot), 

(32) 

pM(x,t)  =  PM(kx,co)exp(ikxx)exp(icot). 

(33) 

Implicit  in  equations  (31)  through  (33)  is  that  the  medium  has  infinite  extent  in  the  x-direction. 
Inserting  equation  (32)  into  equation  (3)  yields 


P(kx,co)  =  Gexp(iyz)  +  exp (-iyz) . 


(34) 


Then,  inserting  equations  (3 1 )  and  (34)  into  equation  (4)  at  z  =  0  gives  the  expression 
P(kx,co)  =  ^-Uz(kx,co)  +  2PI. 


(35) 


The  periodic  mass  term  in  equation  (30)  is  transformed  by  inserting  equation  (31),  applying 
Poisson’s  summation  formula,  and  rewriting  it  as 


D  //  .  _  —  MO)~  V"  TJ  ^ 

p\f(k ;  X  U 


u  —  — oo 


,  2  nn\ 

kx+~— 

L  J 


(36) 


Inserting  equations  (31),  (32),  (35),  and  (36)  into  equation  (29)  yields 


(  2  ^ 
A  T  P  f  G) 

Dkx  -  phor  +— — 


Uz(kx,co)  =  -2PI  + 


Mco 


2  oo 


Z 


n=-oo 


(  2nn\ 

k*+~r 

L  J 


(37) 
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Rearranging  gives 


Uz(kx,co)  =  -2PIY(kx,co)  + 


Mco1 


I  v. 


L/7  =  — oo 


(  2nn\ 
kx+  — 

L  J 


Y(k  ,co), 


(38) 


where 


Y(k  ,co)~ 


4  2  Pf™ 

Dkx  -  phoj  +  -i- 


ly 


(39) 


The  substitution  of 
2nm 


kx  ->  ^  +  • 


(40) 


is  now  applied  to  equation  (38)  yielding 


, r  ,,  2 nm  .  _  _  ....  2nm 

Uz(kx+ - ,co)  =  -2PjY(kx+ - ,©) 

/  / 


+ 


Mo/ 


oo 

,  2nn  2  nm  \ 

I  Uz 

kx+  +  , 

_n=-cc 

v  L  L  J 

....  2^/w 

T(£r  + - ,«). 


(41) 


Taking  the  infinite  summation  of  equation  (41)  with  respect  to  w  gives 


Zr  r  2nm  .  .  n  rn  ....  2nm  . 

Uz(kx+-—,co)  =  -2Pl  X  Y(kx+——,co) 

m=- oo  ^  m=-oo  ^ 


t  Mco2  ^ 


/w=— oo  lL/i=-oo 


oo  ^ 

£  ^ 


2;r«  2;rw? 


*v  + - + 


l  ; 


2nm  .  . 
T(/:r  + - ,o>) 


(42) 


Using  the  shifting  property  of  infinite  summations,  the  second  term  on  the  right-hand  side  of 
equation  (42)  can  be  changed  to 


m=- oo 


00 

( .  Inn  2 Kin  ) 

I  Uz 

**+  + 

- 

_n=-oo 

\  L  L 

y 

00 

II 

M 

k 

W=- 00 

w,  2k  m 
Y{k  + - ,ft>) 


L 

2nn\  00 


I  w, 


2nm 


,co). 


m=- oo 


Inserting  equation  (43)  into  equation  (42)  yields 


Z_  _  2/T/77  __/f  27TA77 

Uz(kx+ - ,co)  =  -2Pj  Y,  Y(kx  + - .®) 

m=—oo  k  m=— oo  k 


Mco 


2  i  i 


n=-oo 


m  -  -co 


Solving  for  the  summated  displacement  term  yields 


-ip,  y  m ,+—,«> 
V  c/z(t, +—,«)  = - a=2 - i - 

£  -  Mco~  2nm 


n=- oo 


1- 


X  K(^+ - ,«) 


/M  =  - oo 


Inserting  equation  (45)  into  equation  (38)  yields 


Uz(kx,co)  =  -2PIY(kx,co)  + 


Mo/ 


-ip,  f.  r(*,+— ,®) 


w=-  00 


.  Mco1  V'  V//  2;rw  \ 

1 - 2-  + - ♦«) 

I  J 

^  m=- oo  L‘ 


After  some  rearrangement,  the  displacement  of  the  system  becomes 


U7(kY,co) 

z  x  =  -2  Y{kxM 


I 


,  Mco~  sr  ,  2/r  m 

L  L 

^  m=- oo  L 


or 
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(43) 


(44) 


(45) 


(46) 


(47) 


>1  -1 


Dkx  + 


A 

Inm 


2 


(48) 


L 


Equation  (48)  is  the  displacement  of  the  system  with  periodic  masses.  Note  that  when  M=  0  or 
L  -»  oo,  the  expression  is  identical  to  equation  (5),  and  when  M  -»  qo,  the  limit  of  the  expression 
is  equal  to  zero. 

Two  different  example  problems  are  now  computed  for  the  thin  plate  response  with  periodic 
masses.  Figure  26  is  an  image  of  thin  plate  normalized  velocity  versus  wavenumber  and  frequency 
generated  using  equations  (28)  and  (48)  for  M=  1 .0  kg/m  and  L  =  0.01  meter.  Figure  27  is  an 
image  of  thin  plate  normalized  velocity  versus  wavenumber  and  frequency  generated  using 
equations  (28)  and  (48)  for  M=  4.0  kg/m  and  L  =  0.01  meter.  Figures  28  through  3 1  are  constant 
frequency  cuts  of  figures  2,  26,  and  27  overlaid  for  a  frequency  range  of  1 .0-4.0  kHz  in  1 .0-kHz 
increments.  (Figures  28  through  31  display  velocity  on  a  decibel  scale.) 
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Figure  26.  Normalized  Velocity  Versus  Wavenumber  and  Frequency  at  M  =  1.0  kg/m 
and  L  =  0.0 1  m — Modeled  Using  Thin  Plate  Theory 
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Figure  27.  Normalized  Velocity  Versus  Wavenumber  and  Frequency  at  M  =  4.0  kg/m 
and  L  =  0.0 1  m — Modeled  Using  Thin  Plate  Theory 
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Figure  28.  Normalized  Velocity  Versus  Wavenumber — Modeled  Using 
Thin  Plate  Theory  with  Masses  at  1.0  kHz 


Figure  29.  Normalized  Velocity  Versus  Wavenumber — Modeled  Using 
Thin  Plate  Theory  with  Masses  at  2.0  kHz 


Figure  30.  Normalized  Velocity  Versus  Wavenumber — Modeled  Using 
Thin  Plate  Theory  with  Masses  at  3.0  kHz 


Figure  3 1.  Normalized  Velocity  Versus  Wavenumber — Modeled  Using 
Thin  Plate  Theory  with  Masses  at  4.0  kHz 
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3.2.  TIMOSHENKO-MINDLIN  PLATE  WITH  MASSES 


The  second  model  investigated  is  that  of  a  thin  plate  with  rotary  inertia  and  shear 
containing  the  periodic  masses  with  a  fluid  load  on  top.  The  equation  of  motion  of  this  system  is 


j4..  ^  ^  ,  3'\  *4 


D 


dx* 


DP.+  Ph~ 


kG  12 


d  uJx,t ) 


2, 3  A  ^4 


dx1dt1 


•  + 


\2kG 


d  u  (x,t)  d  u  (x,t) 
-  +  ph- 


df 


dV 


=  -p(x,t)  + 


D  ' 

P2p(x,t ) 

[ph2] 

VuGh  ) 

K 1 2  kG  j 

PuM. 

dt 


(49) 


where  the  interaction  of  the  fluid  and  the  plate  contains  rotary  inertia  and  shear  effects,  and  the 
interaction  of  the  discrete  masses  with  the  plate  does  not  contain  rotary  inertia  or  shear  effects. 
This  assumption  has  been  included  to  greatly  simplify  the  model.  As  with  the  thin  plate,  the 
displacement  can  be  determined  by 


Pi 


2  Z(kx,G>) 


Mco 2 
L 
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I  Wx 

/w=- oo 


2.nm 

+ - ,co) 
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(50) 


where 
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0>2kl  + 
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12  kG  12iyxG 


p2h3 
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CO 4  + 

( pf  \ 

-P--ph 
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l  J 

(50) 


(O 


Note  that  when  M=  0  or  L  -»  oo,  the  expression  is  identical  to  equation  (9);  when  M  ->  oo,  the 
limit  of  the  expression  is  equal  to  zero. 
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The  two  example  problems  that  were  computed  for  the  thin  plate  response  with  periodic 
masses  are  now  repeated  for  the  thin  plate  with  rotary  inertia  and  shear  effects.  Figure  32  is  an 
image  of  Timoshenko-Mindlin  plate  normalized  velocity  versus  wavenumber  and  frequency 
generated  using  equations  (28)  and  (50)  for  M=  1 .0  kg/m  and  L  -  0.01  meter.  Figure  33  is  an 
image  of  Timoshenko-Mindlin  plate  normalized  velocity  versus  wavenumber  and  frequency 
generated  using  equations  (28)  and  (50)  for  M  =  4.0  kg/m  and  L  =  0.01  meter.  Figures  34 
through  37  are  constant  frequency  cuts  of  figures  3,  32,  and  33  overlaid  for  a  frequency  range  of 
1 .0-4.0  kFlz  in  1 .0-kHz  increments.  (Figures  34  through  37  display  velocity  in  a  decibel  scale.) 
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Normal  Velocity  (dB)  -  Timoshenko-Mindlin  Plate  Model 
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Figure  32.  Normalized  Velocity  Versus  Wavenumber  and  Frequency  at  M  =  1.0  kg/m 
and  L  =  0.01  m — Modeled  Using  Timoshenko-Mindlin  Plate  Theory 
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Figure  33.  Normalized  Velocity  Versus  Wavenumber  and  Frequency  at  M  =  4.0  kg/m 
and  L- 0.01  m — Modeled  Using  Timoshenko-Mindlin  Plate  Theory 
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Frequency  of  1000  Hz 


Figure  34.  Normalized  Velocity  Versus  Wavenumber — Modeled  Using  Timoshenko- 
Mindlin  Plate  Theory  with  Masses  at  1.0  kHz 


Figure  35.  Normalized  Velocity  Versus  Wavenumber — Modeled  Using  Timoshenko- 
Mindlin  Plate  Theory  with  Masses  at  2.0  kHz 
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36.  Normalized  Velocity  Versus  Wavenumber — Modeled  Using  Timoshenko 
Mindlin  Plate  Theory  with  Masses  at  3.0  kHz 


Figure  37.  Normalized  Velocity  Versus  Wavenumber — Modeled  Using  Timoshenko 
Mindlin  Plate  Theory  with  Masses  at  4.0  kHz 


3.3.  THICK  PLATE  WITH  MASSES 


Thick  plate  theory  is  now  investigated  using  the  finite-element  method.  The  problem  of 
added  mass  with  Z.  =  0.0 1  meter  and  M=  4.0  kg/m  is  modeled  using  SARAH2D,  a  commercial 
finite-element  program.14  Initially,  the  plate  is  modeled  without  masses  and  compared  to  the 
results  from  equation  (11).  Figure  38  shows  the  finite-element  analysis  results  for  the  fluid- 
loaded  plate  without  the  masses,  and  figure  39  shows  the  same  results  using  thick  plate  theory. 
Note  that  these  two  figures  are  nearly  identical  except  for  grid  discretization.  Figure  40  is  the 
finite-element  results  for  the  fluid-loaded  thick  plate  with  L  =  0.01  meter  and  M=  4.0  kg/m. 
Figure  41  shows  cuts  of  figures  38  and  40  at  0°  and  30°,  which  correspond  to  an  apparent  wave 
speeds  of  infinity  and  3000  m/s,  respectively.  Notice  that  the  major  effects  of  the  periodic 
masses  are  at  lower  frequencies.  Finally,  figure  42  shows  cuts  of  figures  38  and  40  at  100  and 
1 50  rad/m.  The  masses  have  a  large  effect  on  the  dynamics  of  the  system  at  low  frequency  and 
low  wavenumber — they  have  very  little  impact  on  the  response  of  the  system  at  high  frequency. 
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Figure  38.  Normalized  Velocity  Versus  Wavenumber  and  Frequency 
Without  Mass — Modeled  Using  Finite-Element  Theory 
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Figure  39.  Normalized  Velocity  Versus  Wavenumber  and  Frequency 
Without  Mass — Modeled  Using  Thick  Plate  Theory 
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Figure  40.  Normalized  Velocity  Versus  Wavenumber  and  Frequency  with  M  =  4.0  kg/m 
and  L  =  0.01  m — Modeled  Using  Finite-Element  Theory 
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Figure  41.  Normalized  Velocity  Versus  Frequency — Modeled 
Using  Finite-Element  Theory,  Cuts  at  0°  and  30° 
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Figure  42.  Normalized  Velocity  Versus  Frequency — Modeled 
Using  Finite-Element  Theory,  Cuts  at  100  rad/m  and  150  rad/m 
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4.  ELASTICITY  EQUATIONS  OF  MOTION 


The  derivation  of  the  fully  elastic  equations  of  motion  are  presented  below.  An  earlier 
report7  derived  the  following  stress  relationships  at  the  surface  of  the  plate  for  normal  and 
tangential  stress,  respectively: 


T..(x,b,t)  =  (A  +  2^) 


du 7(x,b,t)  du Jx,b,t) 


dz 


+  A  - 


dx 


=  -px(x,  b,t). 


and 


T  (*,6,0  =  // 


du  (x,b,t)  du7(x,b,t) 


dz 


dx 


=  0. 


The  stress  relationships  at  the  bottom  of  the  plate  are 


T^^^x,u,t)  ( A  +  2//) 


duz(x,a,t)  dux(x,a,t ) 

■  +  A" 


dz 


dx 


=  0, 


and 


Tzx(x,a,t)  =  p 


du  (x,  a,  t)  du7(x,a,t ) 


dz 


=  0. 


(52) 


(53) 


(54) 


(55) 


If  the  mass  is  on  the  bottom  surface  of  the  plate,  equations  (54)  and  (55)  are  modified  to  become 


r..(x,  a,t )  =  (A  +  2/j) 


d 


u.(x,  a,t)  d 
— - +  /t  — 

dz 


ux(x,  a,t) 
dx 


tiAxj) 

dt2 


S(x-nL), 


(56) 


and 


r  (x,  a  =  v 


u  (x, a, t)  d  u,(x,aj) 


dz 


dx 


A  d 
- 

n=-  oo 


dr 


ur(x,t) 

—  S(x-nL).  (57) 
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The  displacements  of  the  system  are 


ux(x,z,t )  =  Ux  (kx ,  z,  a>)exp(ikxx)exp(icot) 


=  yA{k  ,co)ik  exp(/orz)  +  B(k  ,co)ikx  exp (-iaz) 


-  C(kx,co)i/3exp(i/3z)  +  D(k  x,a>)ifl  exp(-i/3z)]exp(ik  xx)exp(icot),  (^8) 


and 


uz{x,z,t )  =  f/z(/:v,z,ry)exp(//::cx)exp(/(uO 


=  [A(k x,co)iaexp{iaz)-  B(kx,co)iaexp(-iaz) 


+  C(kx,co)ikx  exp(i/3z)  +  D(kx,(o)ikx  exp(-/'/?z)]exp(/7:vJc)exp(/Yu/).  (^9) 

To  obtain  a  solution,  equations  (58)  and  (59)  are  inserted  into  equations  (52),  (53),  (56),  and  (57) 
and  the  coefficients  A,  B,C,  and  D  need  to  be  determined.  As  posed,  it  is  unclear  if  this  solution 
exists. 
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5.  CONCLUSIONS 


This  report  has  derived  the  quations  of  motion  of  a  fluid-loaded  thin  plate  that  contains 
embedded  masses  when  it  is  forced  by  an  acoustic  plane  wave  at  a  specific  wavenumber.  The 
equations  were  copied  to  the  fluid  and  the  masses  and  then  manipulated  to  determine  the 
displacement  field.  The  method  was  then  applied  to  a  Timoshenko-Mindlin  plate  and  to  a  thick 
plate. 


Several  conclusions  were  supported  by  this  study.  For  a  soft  elastomeric-type  material 
with  thickness  of  0.04  meter  or  less,  thin  plate  equations  of  motion  are  valid  to  about  1 .5  kHz, 
and  Timoshenko-Mindlin  plate  equations  of  motion  are  valid  to  around  3.5  kHz.  These  results 
are  independent  of  periodic  masses.  The  problem  of  periodic  masses  in  a  thin  plate  and  in  a 
Timoshenko-Mindlin  plate  with  a  fluid  load  on  the  top  and  a  plane  acoustic  wave  excitation  have 
been  derived  and  analyzed  using  different  size  masses.  The  problem  of  periodic  masses  in  a 
thick  plate  with  a  fluid  load  on  the  top  and  a  plane  acoustic  wave  excitation  has  been  analyzed 
using  finite  element  methods  for  one  size  mass.  The  problem  has  also  been  formulated  using 
thick  plate  equations  of  motion. 

It  is  recommended  that  the  thick  shell  equations  of  motion  be  modified  to  handle  the 
periodic  masses  and  examined  to  determine  if  a  closed-form  solution  exists. 
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